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Interrelation between the two-particle and mean-field problems is used to describe the strong pseu- 
dogap and superconducting states in cuprates. We present strong pseudogap state as off-diagonal 
short-range order (ODSRO) originating from quasi-stationary states of the pair of repulsing particles 
with large total momentum (K - pair). Phase transition from the ODSRO state into the off-diagonal 
long-range ordered (ODLRO) superconducting state is associated with Bose-Einstein condensation 
of the K - pairs. A checkerboard spatial order observable in the superconducting state in the 
cuprates is explained by a rise of the K - pair density wave. A competition between the ODSRO 
and ODLRO states leads to the phase diagram typical of the cuprates. Biordered superconducting 
state of coexisting condensates of Cooper pairs with zero momentum and K - pairs explains some 
properties of the cuprates observed below T c : Drude optical conductivity, unconventional isotope 
effect and two-gap quasiparticle spectrum with essentially different energy scales. 
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I. INTRODUCTION 

The high-Xj, cuprate superconductor can be considered 
a doped two-dimensional (2D) antiferromagnetic (AF) in- 
sulator with strong Coulomb repulsion.^ Most commonly, 
one accepts that the physics of the cuprates can be un- 
derstood within the simplest one-band models of strong 
correlated systems such as the Hubbard and t — J mod- 
els. It seems impossible to obtain analytic solutions for 
the ground states of these 2D models. Therefore, a varia- 
tional approach based on a choice of an appropriate trial 
wave function can be considered as a natural way to solve 
the cuprate problem.— 

The wave function of the resonating valence bond 
(RVB) ground stated derived from Gutzwiller projected^ 
d - wave ground state of the Bardeen-Cooper-Schrieffcr 
(BCS) model eliminates a possibility of double occupancy 
of a site (no double- occupancy constraint) and corre- 
sponds to extremely strong on-site correlations. A great 
many of important results, obtained within the RVB ap- 
proach, can be considered as a ground of the physics of 
the cuprates. 

However, the problem is complicated by the fact that, 
in a strong correlated system, ground state energies of 
different ordered states may be close to each others^ d - 
wave superconductor (dSC), staggered flux phase,- spin 
and charge density waves (SDW and CDW, respectively) 
and some others** Account of the competition and co- 
existence of such ordered states within extremely simpli- 
fied models by numerical tools leads to a wide variety of 
phase diagrams. It is clear that, among them, one can 
always find the diagrams reminding those typical of the 
cuprates. 

Various approximations within the on-site Coulomb re- 
pulsion models often lead to antipodal conclusions con- 
cerning a possibility of the superconducting (SC) state 
itselfi 7 - The complete suppression of the double occu- 
pancy under Gutzwiller's projection promotes a rise of 
an insulating rather than SC order. In this connection, 
Laughlin 7 - has proposed an alternative approach to clar- 



ify the problem of the superconductivity of the cuprates. 
Instead of numerical study of a highly simplified Hamil- 
tonian, he suggested to select a reasonable ground state 
in order to determine the Hamiltonian leading to such a 
state. 

To take into account a realistic on-site repulsion, the 
ground state of Laughlin's gossamer superconductor- 
is chosen in the form of an incomplete projected BCS 
d - wave state with partially suppressed double occu- 
pancy. Hamiltonian with such an exact forethought 
ground state, along with strong on-site repulsion, mani- 
fests an attractive term which can lead to dSC stated 

Gossamer superconductivity of the underdoped com- 
pound can be associated with a band of states with rela- 
tively low spectral weight inside a pronounced insulating 
forbidden band so that the chemical potential turns out 
to be pinned near the middle of this band. 7 

The repulsion-induced dSC state is highly sensitive to 
electron dispersion^ The simplest tight-binding approx- 
imation taking into account only the nearest neighbors 
(t - model) seems to be insufficient, therefore, it is nec- 
essary to consider more complicated dispersion with the 
next nearest neighbor terms (t — t' - model). Numerical 
study shows^ that the stable dSC state corresponds to 
hopping integral ratio t'/t within a narrow range near 
t'/t = —0.3. Just the same value of t'/t is consistent 
with available angle resolved photoemission spectroscopy 
(ARPES) data£ 

As follows from the SU(2) approach to the RVB 
problem;? it is necessary to consider doublets of 
fcrmions and bosons to realize the spin-charge separa- 
tion correctly^ Two minima of the SU (2) boson disper- 
sion are relative to the points (0, 0) and (tt, tt) in the 2D 
Brillouin zone^S therefore, the SC pairing channel corre- 
sponding to large pair momentum should be taken into 
account along with the Cooper channel corresponding to 
zero pair momentum. 

Geshkenbein et aljii have assumed that an enhance- 
ment of the scattering between the saddle points of elec- 
tron dispersion results in the fact that the electron- 
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electron interaction with large momentum transfer can 
be "less repulsive" with respect to small transfer. There- 
fore, in the vicinities of the saddle points, fermions may 
be paired into bosons. Such noncoherent preformed pairs 
arising near the antinodal arcs of the Fermi contour 
(FC) might exist in the pseudogap state of underdoped 
cuprates as a normal Bose liquid^ 

The FC outside of the arcs corresponds to unpaired 
fermions coexisting with the preformed pairs. Bosc- 
Einstein condensation (BEC) of the preformed pairs with 
large momentum due to their interaction with unpaired 
particles results in the SC gap on the whole of the FC. 
The SC state that arises in such a way describes reason- 
ably rather wide (intermediate with respect to BCS and 
BEC limiting cases) fluctuation region above T c . 

Instabilities in 2D strong correlated electron system 
were investigated within the t — t' - model at small t' 
by renormalization-group (RG) method o 12 i 13 using a dis- 
cretization of the FC into a finite number of patches. The 
singularity in the Cooper channel exhibits a squared loga- 
rithmic divergence at low energies. For insulating Peierls 
channel with electron- hole pair momentum = (tt, 7r), 
the singularity also exhibits a squared logarithm in the 
particular case t' — when nested FC has the form of 
a square coinciding with the boundary of the magnetic 
Brillouin zone of the parent compound. At t' ^ and low 
doping, that is in the case of a deviation of the FC from 
the perfect nesting, the divergence is found more weak 
with logarithmic enhancement of the order of In \t/t'\ un- 
der the condition that \t'\ <C \t\. The singularities in the 
insulating and SC channels corresponding to zero and 
pair momenta, respectively, are found to be logarithmic 
in the case of small but nonzero \t'/t\. Such a case cor- 
responds to approximately nested FC disposed close to 
saddle-point van Hove singularity. RG approach, involv- 
ing the nesting effects^ gives a possibility to select sin- 
gular contributions into pairing channels but remains cor- 
responding pre-exponential factors to be undetermined. 

General symmetry consideration, based on Zhang's 
SO(5) theory^ or SU(4) theory by Guidry et al.^ shows 
that one should take into account a closed set of compet- 
ing ordered states to describe key features of correlated 
electron system. In this sense, singlet SC pairing channel 
with large momentum incorporating singlet orbital insu- 
lating long-range (possibly, hidden) 17 i 18 or short-range 
(fluctuating between dSC and staggered flux states)^ or- 
der may be naturally connected with the Cooper channel. 
Thus, there should be two SC gap parameters related to 
large and zero pair momentum, respectively. 

The SC gap, which determines T c and corresponds to 
a rise of the coherence in the system of electron pairs, 
can be directly extracted from experiments on Andreev 
reflection 2 ^ or Josephson tunneling^ The observation of 
two SC gaps of about 10 meV and 50 meV, respectively, 
in tunnel experiment 2 ^ in Bi2212 (in particular, a sup- 
pression of the lesser gap in high magnetic field at tem- 
peratures 30 — 50 mK) may be considered as an indi- 
rect evidence in favour of two SC energy scales in the 



cuprates. 

One more energy scale, observed in ARPES and tun- 
nel spectra of underdoped cuprates can be associ- 
ated with the strong pseudogap stated To describe this 
state one can start from a reasonably chosen one-particle 
Green function. Recently, Yang et developed RVB 
phenomenology of the pseudogap state based on the as- 
sumption that this state can be viewed as a liquid formed 
by an array of weakly coupled two-leg Hubbard ladders. 
The coherent part of the Green function obtained within 
the random phase approximation is consistent with the 
Luttinger theorem 2 ^ and describes evolution of the FC 
(from small pockets to closed contour) with doping. Sim- 
ilar results follows from both the spin-charge separation 
approach 27 and phenomenological account of short range 
insulating order above T c £2- 

In this paper, we develop the concept of Coulomb 
pairing 29 leading to the biordered state originating from 
two SC pairing channels, with large and zero pair mo- 
menta. The all-sufficient conditions of repulsion-induced 
superconductivity in these two channels are discussed in 
Sec. II. Sec. Ill deals with the strong pseudogap state 
arising from incoherent quasi-stationary states of pairs 
with large momenta that are inherent in the screened 
Coulomb pairing potential. In Sec. IV, we consider the 
symmetry and two-gap spectrum of the biordered state. 
Finally, some possible manifestations of the strong pseu- 
dogap and biordered SC states are discussed in Sec. V. 



II. COMPETING PAIRING CHANNELS 

Screening of Coulomb repulsion in three-dimensional 
isotropic degenerate electron gas results in momentum 
dependent interaction energy of two electrons, 



U(k) = 47re 2 /[fc 2 e(fc)], 
where static permittivity has the form 3 ^ 
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Here, x = k/2k F , k F and k^ 1 — (4-Ke 2 ng) 1 / 2 are Fermi 
momentum and screening length, respectively, n and g 
are electron concentration and density of states on the 
Fermi level. 

Kohn singularity at k — 2k F leads to the interaction 
energy with damped Friedel oscillation in the real space. 
At a distance r kp 1 , this potential can be written as22, 



U(r) 
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(3) 



Kohn and Luttinger— have argued that attractive contri- 
bution into screened Coulomb repulsion originating from 
Friedel oscillation is sufficient to ensure Cooper pairing 
with non-zero angular momentum. Because of the weak- 
ness of the Kohn singularity, corresponding SC transition 
temperature turns out to be very lowi^I 
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In the case of nested FC, the Kohn singularity trans- 
forms into the Peierls one with strong anisotropy of e(k). 
Therefore, effective pairing interaction can be enhanced 
both in particle-hole and particle-particle channels. In 
particular, this can give rise to CDW or SDW in singlet 
or triplet insulating pairing channels, respectively. 

Peierls singularity in a particle-hole channel originates 
from the fact that momentum transfer turns out to be 
equal to nesting vector Q for any particle on the FC 
(Fig. la). For the sake of simplicity, in Fig. 1 the FC 
is presented as a square corresponding to the t - model 
at half-filling. In the case of low doping and under the 
condition that \t'\ <C t, one can expect that the FC may 
be found close to this squared 

In the Cooper channel, Peierls enhancement of the 
Kohn singularity emerges as appreciably more weak be- 
cause momenta before and after scattering (p and p' 
in Fig. lb, respectively) giving rise to this enhancement 
should be related as p — p 1 « Q. Integration over mo- 
mentum transfer p' — p along the whole of nested FC 
smoothes down the Peierls singularity from the Cooper 
channel. Logarithmic singularity (ln(|£ /e|), where e 
is a cut-off energy) in the Cooper channel is ensured by 
a general feature of electron dispersion, e(— p) = e(p), 
that holds for any momentum p (statistical weight of the 
Cooper pair is proportional to the length of the entire 
FC). 

Density of states of 2D system manifests logarith- 
mic van Hove singularities originating from saddle-point 
vicinities with hyperbolic metric. Due to close proxim- 
ity of the FC and the isoline connecting saddle points 
(±7r, ±7r), effective coupling constant w turns out to be 
logarithmically enhanced, w — > w ■ In (2tk\l\e — £ s |)- 12 
Here, e s is saddle point energy (within the t - model, 
e s = 0) and k c has meaning of a scale of the part of 2D 
Brillouin zone with hyperbolic metric. 

In the case of SC pairing with large total momentum 
K (K - channel), the momenta of the particles compos- 
ing a pair with given momentum (K - pair), both being 
either inside or outside the FC at T = 0, should belong 
to only a part of the Brillouin zone (domain of kinematic 
constraint) rather than the whole one. In a general case, 
kinetic energy of the pair with relative motion momen- 
tum fe, 

2t K (k)=e(K/2 + k) + e(K/2-k)-2fi, (4) 

vanishes only at some points of the FC inside this domain 
(/X is the chemical potential). Therefore, in contrast with 
the Cooper pairing when £g(fe) = on the whole of the 
FC, integration over k eliminates the logarithmic singu- 
larity in the K - channel. However, if kinetic energies, 
e(K/2 + k) and e(K/2 — k) coincide on finite pieces of 
the FC ("pair" Fermi contour, PFC), logarithmic singu- 
larity In (|e /|e|) survives and the K - channel can result 
in the SC order ^ Mirror nesting condition, 
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e{K/2 + k) = e{K/2-k), 
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FIG. 1: Nested Fermi contour (bold line) in the form of 
a square coinciding with 2D magnetic Brillouin zone of 
parent compound with halh-filled conduction band (within 
tight-binding model with nearest-neighbor interactions), a: 
electron-hole pairing with momenta p and p' = p Q (Q 
is nesting momentum); b: Cooper pairing with zero total 
momentum (p and p' are momenta before and after scat- 
tering, respectively); c: SC pairing with nesting momentum 
(p± — Q/2 ± fc, where k is momentum of the relative motion 
of the pair). Domain of kinematic constraint is degenerated 
into a line coinciding with one of the sides of the square; 
d: SC pairing with an incommensurate total momentum K 
(p ± = K /2±fe). Domain of kinematic constraint is bounded 
by the line 1-2-3-4-5-6. Parts 2-6 and 3-5 of this line are 
nested pieces of the FC resulting in singular contribution into 
this pairing channel. 



determines the locus in the momentum space that log- 
arithmically contributes to the K - channel. Statistical 
weight of K - pair, proportional to the length of the PFC, 
can be less in comparison with the Cooper channel. 

Mirror nesting is a necessary (not all-sufficient) con- 
dition of the SC pairing with large momentum. Indeed, 
this condition is perfectly satisfied in the case of nested 
FC when K = (Fig. lc). However, it is obvious that, 
in such a case, domain of kinematic constraint degener- 
ates into a line resulting in zero statistical weight of the 
paired state. To obtain finite statistical weight, one can 
choose incommensurate pair momentum K ^ Q v . This 
results in the domain of kinematic constraint in the form 
of relatively narrow strip containing the PFC as shown in 
Fig. Id. The coefficient of the logarithmic contribution 
to the K - channel should be proportional to the length 
of the PFC. Considering the PFC as two patches con- 
nected by nesting vector Q v (Fig. Id), one can conclude 
that this coefficient has to be logarithmically enhanced 
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by umklapp scattering inherent in the Peierls channel^ 
Patch approximation 12 appears to be relatively good just 
in the case of short PFC. The reason is that integration 
over momentum transfer p — p' w cannot completely 
eliminate enhancement of the pairing interaction due to 
logarithmic singularities of the permittivity. 

Another necessary condition of SC pairing under re- 
pulsion is connected with the existence of oscillating at- 
tractive contribution into the pairing potential. It should 
be noted that an oscillation itself cannot ensure a rise of 
a bound state. For example, simple step-wise repulsive 
potential U(k) = U > defined in a finite domain of 
the momentum space oscillates in the real space. How- 
ever, by analogy with the problem of a bound state in 
one-dimensional asymmetric potential well, 32 such a po- 
tential cannot result in a bound state even under mirror 
nesting. 

One can consider screened Coulomb potential U(k-k') 
as a kernel of Hermite integral operator with complete 
orthonormal system of eigenfunctions defined within do- 
main of kinematic constraint S, 



W 



,(fc) = A a ]T U(k-k')tp 3 (k'). 



(6) 



Here, a set of X s represents the spectrum of a pairing 
operator which can be written in the form of the Hilbert- 
Schmidt expansion, 



A, 



(J) 



The necessary (and sufficient, under mirror nesting) con- 
dition of the SC pairing under repulsion is the existence of 
at least one negative eigenvalue of the pairing operator^ 
In the case of comparatively small domain of kinematic 
constraint, one can replace the screened Coulomb poten- 
tial by its expansion in powers of momentum transfer, 
k = k — k 1 , up to the term of the second order, 



(8) 



where U and r have meaning of an on-site repulsive 
energy and screening length, respectively. The simplest 
repulsive kernel defined inside S, has two even and 
two odd (with respect to inversion k — > — k) eigenfunc- 
tions. Singlet SC order parameter should be determined 
by only even eigenfunctions belonging to eigenvalues X 1 
and A 2 of opposite sign. 

Scattering between nested pieces of the FC leads to 
strong anisotropy of the permittivity. Therefore, expan- 
sion of w(k, k') in powers of momentum transfer close to 
nesting momentum Q appears to be anisotropic as well. 
Resulting pairing interaction kernel, analogous to Eq.®, 
preserves its eigenvalue feature X X X 2 < 0. It should be 
noted that nesting momentum Q ^ therefore, Q has 
the meaning of a new nesting momentum, inherent in the 
real FC, which can result in the Peierls enhancement of 
the SC pairing. 
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FIG. 2: Repulsive pairing potential, w(r), and bound state, 
\ip\ 2 (dotted line), distributions in the real-space (schemati- 
cally). Energies E i and E q correspond to bound and quasi- 
stationary states, respectively. Barrier height E b corresponds 
to a break of the pair without tunnelling through the barrier. 



Matrix of pairing operator Eq. ((7]) between its eigen- 
functions is diagonal, w ss , = Xj 1 5 ss ,. The necessary 
condition of the existence of nontrivial solution to the 
self-consistency equation with kernel ([8]) has the form 
X X X 2 < 0. Written in arbitrary basis, it takes the form 
of the Suhl inequality, 



w n w; 22 - w 12 w 21 < 0, 



(9) 



introduced as a necessary condition of superconductivity 
within two-band model.— 

The effective pairing potential oscillates in the real 
space and, in agreement with Laughlin's proposal^ man- 
ifests repulsive core at small distance (corresponding to 
incomplete no double-occupancy constraint) and attrac- 
tive contribution outside of the core. Thus, there is a 
possibility of a rise not only of a bound (with negative 
energy, E < 0) but also of quasi-stationary (with E > 0) 
paired state with large momentum (Fig. 2). 

The singular contribution into SC order parameter is 
determined by relatively small vicinity of the PFC with 
an energy scale e . In this respect, repulsion- induced K - 
pairing seems to be similar to phonon-mediated pairing 
arising from attraction with negative coupling constant 
V. Rough estimation^! of Coulomb repulsion within the 
phonon-mediated mechanism of superconductivity leads 
to the fact that, to ensure SC pairing, |V| must exceed a 
threshold value, 



\V\ > 



l+gUJn(E F /e D )- 



(10) 



Here, E F and U c are Fermi and average Coulomb ener- 
gies, respectively. Phonon-mediated attraction is defined 
inside a narrow layer (domain of dynamic constraint with 
energy scale of the order of Debye energy e D ) enveloping 
the FC. 

One can expect that a deviation from perfect mirror 
nesting condition ((5|) outside of the PFC does not elimi- 
nate the logarithmic singularity in the K - channel. Typ- 
ical of high - T c cuprates, nearly nested FC in the form 
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of a square with rounded corners^ and finite screening 
length of the Coulomb interaction result in a weak sup- 
pression of the singularity. The SC gap turns out to be 
preserved on the PFC (by analogy with the case of insu- 
lating pairing studied by Losovik and Yudson)^ 

The Cooper and K - channels can be found as approx- 
imately equally enhanced by both Peierls singularity of 
screening and proximity of the FC to van Hove singular- 
ity of the density of states. Therefore, it is the kinematic 
constraint that can ensure the preference of the K - chan- 
nel over the Cooper channel. One can believe that, in 
doped cuprates, the singular contribution to Peierls en- 
hanced Coulomb K - pairing can dominate the Cooper 
channel. 



III. STRONG PSEUDOGAP STATE 

Taking into account the ground state instability due to 
a rise of pairs, the mean-field approach to the problem of 
superconductivity excludes fluctuations of paired states 
from consideration. Within the mean-field theory^ the 
SC gap is directly relative to the binding energy of a 
pair resulted from the two-particle Cooper problem.— In 
the case of the K - pairing, such a problem may admit 
more complicated solution as compared to the attraction- 
induced pairing.— 

Integral equation which determines a wave function of 
the relative motion of two interacting particles (holes) 
above (below) the FC can be written as 



ijj(k) = G (0) (w;fc) w(k,k')^(k') 



(11) 



fc'GH 



Here, lo and k (k') are an energy and momentum of 
the relative motion before (after) scattering, respectively, 
and 



G (0 V;fc) = [w-C(fc)+*7]" 



(12) 



is one-particle Green function corresponding to free rel- 
ative motion of the K - pair, 7 — > +0. In contrast 
with one-particle Landau Fermi liquid Green function, 
the condition that [G^ (0, fc)] _1 = does not determine 
a closed FC. Indeed, a locus in the momentum space 
resulting from this condition written in equivalent form 
£ K (k) — 0, is either some isolated points or finite pieces 
of mirror nested FC. 

In the case of mirror nested FC, one can separate a 
singular contribution to the Green function originating 
from relatively small (with energy scale e ) part S s of 
the domain of kinematic constraint 5. The rest of the 
domain, including an energy range from e up to a cut-off 
value of about E F , results in a regular contribution into 
G<°). One can consider this contribution in a way similar 
to the account for the Coulomb repulsion within phonon- 
mediated pairing attraction scenario^ and renormalizc 
kernel Eq. to a kernel, defined inside S s , with the 
same spectrum. 



One-particle Green function G(lo; k) corresponding to 
relative motion of K - pair of particles (holes) excited 
above (below) the FC can be represented in the basis 
formed by the eigenfunctions of the renormalized pairing 
operator w(k, k'), 

G ss ,M = ^:(fc)G(^;fcV s ,(fc). (13) 
fees 

Matrix elements (fT5|) are the solutions to Dyson equation, 
J2{5 SS „ \-,}gW(u)}G s ,,A") = GfjM, (14) 

s" 

in which matrix elements G^l (lo) of free Green function 
(fT2|) are defined similar to Eq. tfH]) . 

Pairing operator with two even eigenfunctions^, results 
in 2 x 2 matrix (TT5)) . One can resolve Eq. ([Tl| with respect 
to G ss , (lo) and then obtain G(ui; k) in the form 



G(w;k) = D- 1 (w)[G^{u;k)-B(w;k)], (15) 



where 



B(u>; k) = Ar^BM J2 Kl^sM 2 , (16) 

s=l 

B(w) = G^MG^M - G^HG^M, (17) 



and 



G^M G£'M , B(w) 



(0), 



D(lu) = 1 - 



A, 



AiA 2 



(18) 



In the case of mirror nested FC, Green function (jl5|) 
manifests a pole resulting in a bound state with negative 
energy lu — E i determined from equation D(lj) = in 

which all functions G^l (u>) are real. This pole is related 
to instability of the ground state with respect to a rise of 
pairs. Within a small vicinity of the pole, Green function 
(fT5)l can be represented as 



[G^(u;;k)-B(u;k)] 1 

G(w;fc) = — im — ^ 



(19) 



where D' = d D/d lo. 

At to > 0, Green functions (|13|) are complex and equa- 



tion D(lo) = can lead to complex solution lo = E q — iT 
where E and Y have meanings of energy and decay of 
quasi-stationary state (QSS) of the relative motion of K - 
pair, respectively— Near this complex pole, Green func- 
tion (TT5)) has the form of Eq. (|19[) where E i should be 
replaced by E q — iY. 

Wave functions of the relative motion of K - pair cor- 
responding to both bound state and QSS, are localized, 
in main, in a wide region of the real space outside the 
repulsive core as shown schematically in Fig. 2. 
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The K - pairs can exist above T c as long-living QSS 
due to considerable increase of density of states in a nar- 
row vicinity of E q . To overcome the potential barrier 
before tunnel decay, such a non-coherent pair should ac- 
cumulate an energy exceeding barrier height E b . Thus, 
the energy E q — E i is sufficient to destroy SC coherence 
whereas corresponding pair-break energy should exceed 
E b — E { . A temperature range between the SC tran- 
sition temperature T c ~ E q — E i and a crossover one, 
T* tr ~ E b — E i , can be interpreted as a strong pseudo- 
gap state observable above T c in underdoped cuprates. 
If density-of-states peak at u> w E q turns out to be 
smoothed due to T being large enough, the strong pseu- 
dogap state becomes unobservable. In such a case, the SC 
transition from coherent into non-coherent state should 
be assisted with a break of pairs at energies « E b — E i 
similar to that in the BCS theory. 

By analogy with the interrelation between the Cooper 
two-particle problem 3 -! and BCS theory^ the pair-break 
energy E b — E i due to direct excitation of particles from 
bound state into continuous spectrum should be trans- 
formed into momentum-dependent energy gap A(fe) in 
the quasiparticle spectrum. In the strong pseudogap 
state, this gap, due to a non-coherence of QSS, can be 



presented as A = yj A;? + A 2 . Here, A c ~ E q — E t corre- 
sponds to transition from the coherent into non-coherent 
QSS and A p ~ E b — E q can be related to a break of If - 
pair as a result of transition between two non-coherent 
states. 

Microscopically, SC gap A c and strong pseudogap A p 
emerge with random phases. Therefore, mean-field value 
A p vanishes at any temperature whereas A c becomes 
nonzero below T c due to Bose condensation of K - pairs 
from QSS into the bound state and vanishes only above 
T c . However, nonzero mean square strong pseudogap, 
|A p | 2 7^ 0, may become apparent well above T c . In this 
sense, pseudogap parameter A p , corresponding to decay 
of QSS of K - pairs, reminds RVB spin liquid pseudo- 
gap introduced by Yang et al. 25 However, it has different 
physical meaning. 

Green function G^ (0; k) changes sign on the PFC 
from positive to negative through an infinity. Therefore, 
Green function (fT5|) at u> = manifests the same feature. 
It should be noted that, in the case of Cooper pairing, the 
Green function changes sign on the whole of the FCj2£ 
In addition, Green function G^ (0; fcj changes sign on a 
zero line determined by equation G'°^(0;fc) = B(0;k). 
This line does not coincide with the PFC. 



Green function (|15[) of the two-particle problem has a 
pole corresponding to a bound state of the relative mo- 
tion of the pair. Therefore, one can suppose, in line with 
Yang et al.f 2 ^ a phenomenological BCS-like form of the 
coherent contribution to the normal (diagonal) Gor'kov 
Green function of the mean-field problem^ 9 . 



G{u;k) = z k 



_{k) 



-(*) 



E{k) + iY 



E{k) - iT 



where E = + |A| 2 and 2u± — 1 ± i K jE are quasi- 
particle energy and coherence factors, respectively. In ac- 
cordance with [fT5|) . one should suppose that momentum- 
dependent quasiparticle weight z k vanishes on the line of 
zeroes and corresponds to a finite value z (0 < z < 1) on 
the PFC. Two terms can be referred to pairs above and 
below the FC. 

Diagonal Green function ([2H|) describes a non- 
superconducting state with off-diagonal short-range 
(ODSRO) corresponding to the existence of non-coherent 
K - pairs above T c . Below T c , the ODSRO transforms 
into the off-diagonal long-range order (ODLRO) intro- 
duced by Yang 3 

Excitation with a transition from the bound paired 
state into long-living QSS corresponds to quite small but 
finite decay T = T(uj;k). The transitions into station- 
ary states above barrier energy E b should be associated 
with an infinitesimal decay, 7 — > +0, leading to conven- 
tional Fermi-liquid behavior of diagonal Gor'kov func- 
tion d^UJ) above T* tr . Thus, a rise of QSS results in 
a non-Fermi-liquid behavior of diagonal Green function 
(|15p that can be manifested in rather wide temperature 
range T c < T < T* tr relating to strong pseudogap state. 
This range corresponds to transitions between boson-like 
bound and quasi-stationary states. Therefore, Eq. (|2"0|) 
can be considered as a bridge between the BCS and BEC 
approaches to the problem of superconductivity, in accor- 
dance with the assumption by Geshkenbein et alM- 



IV. SUPERCONDUCTING STATE WITH 
LARGE MOMENTUM 

One can believe that the nearly nested FC of under- 
doped (up to optimum doping) cuprate compound in the 
form of a square with rounded corners, shown schemat- 
ically in Fig. 3, results in the fact that the K - chan- 
nel corresponding to an incommensurate momentum K 
dominates the Cooper channel. In such a case, it is a 
rise of coherence in the system of K - pairs that deter- 
mines SC transition temperature T c , Therefore, there is 
a temperature range well below T c in which mean-field 
SC order parameter A c (fe), relating to the K - channel, 
can be approximately considered as governed by the only 
self-consistency equation, 



A c (fc) 



iM-fBl. (2i) 



Quasiparticle energy, 



E(k) 



e K (k) 



|A c (fc)|* + |A p (fc)P 



(22) 



(20) 



aside from A c (fe), includes strong pseudogap parame- 
ter, A p (fc), associated with QSS that can exist above 
T c . Thus, Eq. (f2Tj) reflects the fact that SC order arises 
from the state other than normal Fermi liquid. There- 
fore, Eq. (f2"Tj) differs from the conventional BCS self- 
consistency equation. It is reasonable to assume that 
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FIG. 3: Schematic representation of the Fermi contour of un- 
derdoped cuprate superconductor, in accordance with Ref.— 
Nested (with nesting momentum Q) and mirror nested (cor- 
responding to total pair momentum K) pieces of the FC are 
shown by bold lines. Dotted line shows magnetic Brillouin 
zone boundary. Shaded: dark narrow ovals correspond to 
vicinities of Peierls singularity of screening; light half-ovals 
designate extended saddle point vicinities. 



A p (fe) is small above optimum doping and gradually in- 
creases with underdoping. This leads to smoothing of 
the singularity in Eq. (f2"Tj) and, as a result, to a gradual 
decrease in T c with underdoping. 

One can reduce summation in Eq. (|2ip to small part 
S s of the domain of kinematic constraint similar to the 
two-particle problem considered above. Under repulsive 
interaction, a non-trivial solution to Eq. (|2"Tj) can arise 
due to a competition of positive and negative contribu- 
tions to the right-hand side of this equation. Thus, such 
a solution should have a line of zeroes (nodal line, NL) 
inside S s . 

In the case of phonon-mediated pairing with account 
of Coulomb repulsion^ the NL of the order parameter 
coincides with the boundary enclosing the domain of dy- 
namic constraint. In this domain, attraction dominates 
logarithmically weakened repulsion in accordance with 
Eq. (fTOjl . This NL is disposed everywhere outside of the 
FC, therefore, mirror nesting of the FC can be considered 
as the only condition of the K - pairing under attraction. 

Peierls enhancement of the pairing interaction results 
in a strong anisotropy of the NL disposed close to the 
FC inside 3 S . This corresponds to an increase in the 
dominant part of S s that mainly contributes into the 
logarithmic singularity in the self-consistency equation 
and thus increases a magnitude of A c (it is clear, that 
|A C | should be much lesser than e ). 

Since the SC order parameter arising in the K - chan- 
nel is essentially momentum-dependent, there are three 
characteristic lines of zeroes: 1) the PFC on which ki- 
netic energy of the pair equals zero, 2£ R -(fc) = 0; 2) the 
NL of the order parameter determined by A c (fc) = 0; 3) 



the curve on which quasiparticle group velocity changes 
sign, V k E(k) = 0. 

These three lines may have common points of intersec- 
tion inside S s , therefore, the NL can be disposed both 
above and below the PFC. This results in qualitatively 
different non-monotonic momentum dependence of co- 
herence factors u±(fe) for two kinds of directions in the 
k - space intersecting at first the PFC and then the NL 
and vice versa. Under pairing repulsion, the scattering 
across the NL turns out to be dominating in comparison 
with scattering inside or outside the NL in accordance 
with Suhl inequality Eq.©. 

Due to the fact that |A p | 2 ^ in the strong pseudogap 
state, coherence factors in diagonal Gor'kov function PO)) 
may overlap each other near the PFC even above T c . On 
the contrary, the BCS coherence factors are step-wise 
functions without an overlap in the normal Fermi liquid 
state. 

The SC state that arises below T c should be described 
by both diagonal and off-diagonal (anomalous) Gor'kov 
functions. Taking into account the fact that mean- field 
(averaged over random phases) pseudogap parameter A p 
vanishes whereas mean-field SC condensate parameter 
A c 7^ below T c , one can introduce off-diagonal Gor'kov 
function F + (uj;k) in a phenomenological way similar to 
that we use to obtain diagonal Gor'kov function (|20[) . 
This function describes the ODLRO stated and can be 
written as 



F + (uj;k) = - 



(lu - E{k) + iT)(u + E{k) - iT) ' 



(23) 



Factor z k is defined inside each of the crystal equivalent 
domains of kinematic constraint EL where j — 1,2,3,4 
in the case of tetragonal symmetry of cuprate planes. 
Paired states with large total momenta Kj , both coher- 
ent and non-coherent, arise exactly inside these domains. 

Parameters A„„-, A_- and r„- are identical for any of S„ 

c j 1 pj j _ J j 

differing only by the domain of definition of k. 

In the whole of the Brillouin zone, the SC order pa- 
rameter in the mixed representation, A ■(!?, fe), can be 
presented as a superposition, 



A c (R,k) 



4 



7 ,(fc)e iK * R A c ,(fc), 



(24) 



where R is center-of-mass radius- vector, coefficients 
Tj-(fe) should be chosen in accordance with the symmetry 
of the order parameter. 

As a function of R, order parameter (|24p arising in the 
K - channel turns out to be spatially modulated similar 
to Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state. 42 i 43 
Such a modulation with relatively short wavelength re- 
flects a variation of K - pair density and can be asso- 
ciated with a pair density wave with a checkerboard or- 
der in cuprate planes. Thus, the K - pairing leads to 
a microscopic ground of the pair density wave concept 
introduced phenomenologically by Zhangi 1 ^ 
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The SC order parameter can be expanded over the 
eigenfunctions of pairing operator (J7J), 29 

A c (R,k)=J2 A cs(R)<Ps(k). (25) 

s 

In the simplest case, repulsion-induced SC pairing can 
be described by two-component order parameter, so that 
two complex components, A cs (i2), s = 1, 2, form the or- 
der parameter structure in the framework of Ginzburg- 
Landau phenomenology. Absolute values and relative 
phase of the components are connected with the relative 
motion of K - pair. 

As follows from Ginzburg-Landau equation system for 
two-component order parameter^ two qualitatively dif- 
ferent SC states become admissible. One of them corre- 
sponds to constant value it of the relative phase of the 
components that is generic for repulsion-induced super- 
conductivity. The other state, with the relative phase 
different from n, can be related to a change of the phase 
of the wave function of K - pair due to a rise of inter- 
nal magnetic field of spontaneous orbital currents. Such 
currents can be associated^ with insulating orbital anti- 
ferromagnetic order, for example, in the form of DDW. 18 



V. BIORDERED SUPERCONDUCTING STATE 

Superposition (j24|) mixes, in particular, two paired 
states with opposite momenta, K and — K. It is clear 
that particles composing pairs with these momenta can 
compose pairs with zero momentum. Thus, the Cooper 
channel appears to be associated with the K - channel 
in accordance with symmetry consideration. 

In such a case, mean- field order parameters A (fc) and 
A c (fe) corresponding to the Cooper and K - channels, 
respectively, should be the solution to a self-consistency 
equation system. This system degenerates into two in- 
dependent equations (each of them determines one of 
the order parameters) if one neglects the interconnec- 
tion of the channels. Then, one would obtain temper- 
atures, T c and T' c , of transitions into the states with 
order parameters A c (fe) and A (fe), respectively. Let 
us assume that T' c < T c even in the case when attrac- 
tive phonon-mediated pairing contributes to the Cooper 
channel. Then, SC transition temperature T c can be ob- 
tained directly from Eq. (|2ip . 

At T' c < T < T c , there arises SC order due to K - 
pairing with order parameter A c (fe) defined in relatively 
small vicinities of the PFC where factor z k is close to 
unity. Inside this temperature range, Cooper pairing on 
the whole of the FC should be induced by K - pairing. 
In this case, the magnitude of A (fc) has to be small in 
comparison with the magnitude of A c (fe). 

Thus, biordered SC state arising in such a way should 
be described by two order parameters, A c and A , de- 
fined in the vicinities of the PFC and entire FC, re- 
spectively. As temperature decreases from T c down to 



T w T' c , Cooper ordering with order parameter A ex- 
ists as induced by the K - channel of SC pairing. In this 
case, the superfluid density turns out to be approximately 
proportional to the PFC length. Opening of the Cooper 
channel at T ss T' c leads to a considerable increase in A 
and, as a result, in the superfluid density which becomes 
proportional to the whole of the FC length at T <T' C . 

In the vicinities of the PFC, two branches (m = 1, 2) of 
strong anisotropic quasiprticle spectrum of the biordered 
superconductor can be written in the form 

E m (k) = ^(fe) + |A p (fe)| 2 + |A c (fe)±A (fc)p. (26) 

Here, we take into account the fact that kinetic energy 
of K - pair (j3J) is equal to kinetic energy of Cooper pair, 

2£ (fc) = e(K/2 + k)+ e(-K/2 - k) - 2p. (27) 

Two-gap spectrum (f26|) . with the lesser gap |A C — A | 
observable at excitation energies up to the greater gap 
A c + A |, should be apparent at T < T' c . Above the 
greater gap the spectral weight transfers from the low to 
high-energy branch of quasiparticle spectrum. 

Diagonal and off-diagonal Gor'kov functions of the 
biordered SC state preserve their form, Eq. (|2"0]) and (|23[) . 
respectively, with the exception of the fact that SC order 
parameter A c has to take into account both SC pairing 
channels. Thus, these two channels result in two coexist- 
ing ODLRO states. 

Most likely, Cooper channel, including both Coulomb 
and phonon-mediated pairing, cannot result in a rise of 
QSS. Therefore, SC gap parameter turns out to be BCS- 
like everywhere on the FC with the exception of the PFC 
on which the unconventional K - channel is opened. 

Symmetry of the biordered SC state is determined by 
Eqs. ([24|) and (|25|) where A s should be considered as the 
components of momentum dependent order parameter 
arising as a result of both Cooper and K - pairing. One 
can approximately represent the gap parameter in the 
conventional form 

A(k) = D(k)(cos k x ± cos k y ) (28) 

where momentum dependent magnitude D(k) reflects 
mainly strongly anisotropic contribution of the K - pair- 
ing. Upper (lower) sign in Eq.([28|) corresponds to ex- 
tended s(d) - wave symmetry of the order parameter. 
It should be noted that, in the case of Peierls enhanced 
K - pairing, nodal line of order parameter A(fe) may 
pass through the center of the corresponding domain of 
kinematic constraint. This results in the fact that order 
parameter A(fe) differs in sign on the opposite parts of 
the PFC. 



VI. CONCLUSION 

We believe that biordered superconductivity may be 
generic for such superconductors as doped cuprate com- 
pounds. Unconventional features of these compounds, 
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FIG. 4: Schematic phase diagram of underdoped cuprates. 
Here, T N and T c are phase transition temperatures corre- 
sponding to Neel AF and unconventional SC orders, respec- 
tively; Tstr is strong pseudogap (sPG) crossover temperature, 
T* bounds the region of weak pseudogap (wPG) state. High 
doping and temperature state corresponds to normal Fermi 
liquid (FL). On top, there are shown doping ranges in which 
one can expect a rise of QSS, QSS together with BS and BS. 



especially, universality of their phase diagram (Fig. 4), 
can be associated with evolution of the FC and pairing 
interaction with doping. 

It is clear that singular contribution to Coulomb pair- 
ing interaction is sensitive to doping dependent form of 
the FC. Therefore, one can suppose that oscillating real- 
space pairing potential w(r) varies with doping in such a 
way that only noncoherent QSS of SC pairs arise under 
extremely low doping. This corresponds to strong pseu- 
dogap penetrating into insulating region of doping below 
the onset of superconductivity at x = x^. In underdoped 
region < x < x opt , along with QSS, there is a bound 
state. Both bound state energy \Ej\ and QSS decay T 
increase with doping so that, near optimal doping x opt , 
pair-break energy approximately coincides with the en- 
ergy corresponding to the loss of phase coherence. Thus, 
in overdoped regime x opt < x < x* , pairing interaction 
can result in only bound state. 

As doping increases, the K - channel may be domi- 
nated by the Cooper channel and overdoped SC state can 
manifest properties inherent in conventional BCS state. 
When doping exceeds x opt , a decrease in T c down to 
T c = at x — x* can be also associated with doping 
dependence of the pairing interaction. This interaction 
becomes more repulsive at x > x opt as the FC gradually 
leaves the vicinity of the extended van Hove saddle point. 

A suppression of the phonon-mediated component of 
the SC pairing may have the same origin, Effective in- 
crease in repulsion can result in the fact that inequality 
(fTU)) can be reversed because of not too large ratio E F /e D 
typical of cuprates. 

There is rather strong evidence that, in underdoped 
cuprates, dimensionless ratio 2A^°^/T C considerably ex- 



ceeds universal BCS value 3.524S. Here, A^ ) has mean- 
ing of the SC energy gap extrapolated down to T = 0. 
In underdoped biordered superconductor, this parameter 
should be determined by both Cooper and K - pairing, 

therefore, one can assume that A^ ' = y^A 2 + A 2 , + A 2 . 

Taking into account the fact that T c is determined by 
K - pairing only and A 2 , A 2 , and A 2 are, generally 
speaking, of the same order, one can easily conclude that 
ratio 2A(°'/T C may considerably exceed 3.52 (observed 
in Ref4^ values 2A^ /T c > 10). In overdoped regime, 
strong pseudogap parameter A p — * and the Cooper 

channel dominates K - pairing. Therefore, 2A<°) /T c 
should be close to 3.52 in accordance with the BCS the- 
ory. 

Superfluid density p s should be determined by con- 
densation of K - pairs within a broad temperature range 
below T c , down to the onset of the Cooper channel. Be- 
low T' c , superfluid density increases considerably; con- 
versely, off-condensate particle density decreases. Drude- 
like behavior of the coherent contribution into optical 
conductivity <j 1 (oj) ~ w~ 2 , observed below T c , can be 
connected with rather high off-condensate density (ex- 
perimental data available^ show that spectral weight of 
the off-condensate particles may exceed spectral weight 
of the SC condensate below T c ). At T < T' c , Drude 
component of the optical conductivity of biordered su- 
perconductor, <Ji(u>), should be suppressed due to shed- 
ding of the off-condensate particles into the condensate 
of Cooper pairs. 

Two-gap excitation spectrum of biordered supercon- 
ductor should be consistent with tunnel conductance 
measurements^ and quasi-linear temperature depen- 
dence of heat capacity, c v = 7(T)T4£ 

All homologous cuprate series investigated demon- 
strate universal dependence of T c on the number of 
Cu02 layers in the unitary cell, T c (n), with maximum 
at n = 34^ Strong initial increase in T c (n) cannot be as- 
sociated with local real-space pairing interaction. Weak 
interlayer tunnelling can explain this feature qualita- 
tively by rather small effective enhancement of the cou- 
pling constant^ Coulomb pairing with finite screening 
length ensures strong correlation between electrons in the 
nearest-neighbor layers and results in a quantitative ex- 
planation of T c (n) leading to almost triple increase of the 
coupling constant^ 

Doping dependent isotope effect, observed in cuprate 
superconductors^ is highly sensitive to sample quality 
and reflects the contribution of phonon-mediated com- 
ponent to the SC pairing interaction. Depending on 
the interrelation between Coulomb and phonon-mediated 
contributions^ the exponent of the isotope effect on T c 
can be close to both zero, in the case of dominating repul- 
sion, and BCS limit in the opposite case of dominating 
phonon-mediated attraction. Relative isotope shift, neg- 
ligible above optimum doping, increases with underdop- 
ing. We believe that this can be considered as an indi- 
rect evidence in behalf of the fact that, in the case of low 
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doping, Coulomb correlation effects dominate phonon- 
mediated contribution to the K - channel which deter- 
mines T c in biordered superconductor. 

The isotope effect on the London penetration depth 
X L , absent within the BCS theory, also turns out to be 
enhanced with underdopingj— The penetration length is 
weakly sensitive to isotope substitution in a wide temper- 
ature range well below T c . Then, starting from T w T' c , 
the isotope shift on X L increases gradually at T — > 0. As 
AT, 2 ~ p s , such a behavior of isotope effect on \ L can be 
associated with temperature and doping dependence od 



superfluid density inherent in the biordered SC state. 
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